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Status report is made of our quenched study of heavy-light matrix elements employing the Wilson quark action 
for heavy quark. Results obtained up to now with 200 configurations at (3 — 6.1 on a 24 3 x 64 lattice and with 
100 configurations at /3 = 6.3 on a 32 3 x 80 lattice suggest that the pseudoscalar decay constant varies little over 
this range of j3 in both charm and bottom regions. Results for the B parameter are also reported. 



1. Introduction 

The B meson decay constant is an experimen- 
tally important quantity that can be calculated 
within the present techniques of lattice QCD. The 
calculations, however, would be hampered by a 
number of systematic effects and problems, which 
do not appear in the calculation for light quarks. 
We have been carrying out a set of quenched sim- 
ulations covering the range of inverse lattice spac- 
ing a -1 = 2 — 4 GeV employing the Wilson quark 
action to investigate systematic errors in the cal- 
culation. In particular, our prime aim is to see 
the effect of scaling violations in the decay con- 
stant associated with heavy quarks of 0(mQa), 
which allows various definitions of meson mass, 
and to examine whether the results continue to 
the static value at the infinite quark mass limit. 
In this article we report results obtained so far at 
= 6.1 and 6.3 using VPP500/80 at KEK. 

2. Run Parameters and fitting procedure 

In Table | we summarize parameters of our 
runs. The analyses reported in this article is 
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based on data sets, the statistics of which are dou- 
bled since the time of the Symposium. 

We choose the spatial lattice size so that the 
physical lattice size is kept approximately con- 
stant at La w 1.9fm, which should be reason- 
ably large to avoid finite-size effects for analyses 
of heavy-light meson properties. The temporal 
lattice size are taken large in order to ensure a 
reliable extraction of the ground state. Gauge 
configurations are generated with the 5-hit pseu- 
doheatbath algorithm separated by 2000 sweeps 
for (3 = 6.1 and 5000 for (3 = 6.3. 

We take seven values of the heavy quark hop- 
ping parameter Kq to cover the mass of charm 
and bottom quarks. Four values of light quark 
hopping parameter K q are employed for an ex- 
trapolation to the chiral limit. In the heavy- light 
sector we calculate meson propagators for an ar- 
bitrary combination of K q and Kq. Light hadron 
propagators are evaluated only for the degenerate 
case. 

Quark propagators are obtained with the 
red/black minimal residual algorithm for the 
point and wall sources in the Coulomb gauge, and 
local operators are used for sink. 

For x 2 fitting of hadron propagators we employ 
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Table 1 

Parameters of runs. 





6.1 


6.3 


size 


24 3 x 64 


32 3 x 80 


K q 


0.1520-0.1543 


0.1500- 0.1513 


Kq 


0.070- 0.142 


0.086- 0.140 


#conf. 


200 


100 


L (fm) 


1.88(7) 


1.89(8) 


a-^GeV) 


2.56(9) 


3.38(15) 


/*(MeV) 


122(3) 


131(4) 


A<°> (MeV) 239(8) 


246(11) 


K c 


0.15499(1) 


0.15183(1) 


strange 


0.1531(1) 


0.1505(1) 



a single exponential uncorrelated fit with a single 
elimination jackknife analysis to estimate errors 
of masses and amplitudes. For the fitting interval 
[t m in.,tmax], we take [14,26] at (3=6.1 and [16,32] 
at (3 = 6.3. These are chosen by the conditions 
that the point and wall source results for effective 
mass overlap at t m i ni and that t max should be as 
large as possible until errors become unacceptably 
large and the fit becomes unstable. 

3. Scale setting 

We employ the p meson mass extrapolated to 
the chiral limit to fix the physical scale of lattice 
spacing. Results obtained with the wall source 
are listed in Table [l] together with the critical 
hopping parameter K c and the value of K corre- 
sponding to the strange quark determined from 
mif/m p . Values for the point source are con- 
sistent, albeit with an error of order 50% larger 
for a -1 . The lattice spacing at (3 = 6.1 is 
consistent with previous estimates in the range 
p re 6.0 - 6.20], while that at = 6.3 is 10% 
higher than the value reported in Ref. ||. 

We also list in Table [l] our results for the pion 
decay constant /„• using the physical scale as es- 
timated above. For the renormalization factor we 
employ the tadpole improved one-loop result || 
Z A =(1- 0.31a v (l/a))/8K c which equals 0.750 
at 0=6.1 and 0.772 at £=6.3. We observe that 
the value of /„• is two standard deviations smaller 
than the physical value at (3=6.1 in contrast with 
the case of /3=6.3 where they are consistent. 
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Figure 1. Kinetic mass as a function of pole mass 
at (3 = 6.3 for heavy-light meson. Light quark 
mass is fixed at K q = 0.1505. Solid curve repre- 
sents the relation TOkin = sinhTO po i e . 

Finally the QCD A parameter A^L estimated 
from the tadpole-improved A/5 coupling constant 
indicates the presence of a small violation of 
asymptotic scaling in this quantity, as has been 
known from previous workjlJ. 

4. Heavy meson kinetic mass 

A manifestation of O(mga) effects for the Wil- 
son quark action appears in the discrepancy of 
the pole mass to po 1c and the kinetic mass m^m 
defined by an expansion of energy of the form, 



E = m 



pole 



P 



P 



2m kil 



' 8 ™kin 



(1) 



where p 2 = J^ 3 =1 (2 sinpi/2) 2 and the constant c 
parametrizes relativistic corrections. It has been 
argued |jj that dynamical scales of heavy-light 
mesons are controlled by the kinetic mass, and 
hence it is a physically more relevant mass pa- 
rameter than the pole mass. We adopt this view 
and examine the mass dependence of the decay 
constant in terms of the kinetic mass. 

In Fig. we illustrate our data for the two 
masses. The kinetic mass is extracted by a fit of 
the measured dispersion relation using four val- 
ues of momenta up to (1,1,1) in terms of (|l|). It 
is interesting to observe that the relation of the 
two masses are reasonably described by the form 
Wkin = sinhmpoie which follows from the disper- 
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Figure 2. Ratio of correlation functions in (g) 
for point (open circles) and wall (filled circles) 
source for K Q =0.086 and # 9 =0.1505 at (3=6.3. 
Solid and dashed lines are fit value and error as 
explained in text. 

sion relation for a free scalar particle given by 
(2sinh(-j)) =mo+p 2 [|). The relation, however, 
is not accurate in the region of charmed meson. 

5. Heavy-light decay constant 

Our method to obtain the lattice value of the 
heavy-light decay constant is summarized by 

(fpV^) lat = V$ K c - 6K q ^8K c - 6K Q 

xZ A ? p V2ZP. (2) 

The first two factors represent the non-relativistic 
normalization for the quark wave functional). 
The constant Z A l p is extracted from the relation, 



(0|A 4 (t)P(0)t|0) t »i 



(0|P(t)P(0)I|0) 



Z A/p tanhm po i c (t- J L f /2).(3) 



Here A4 and P at time t are the local axial-vector 
and pseudoscalar currents, the subscript s at time 
t = denotes the type of source, and the hy- 
perbolic tangent takes into account the change of 
sign of the numerator at half the temporal lat- 
tice size Lt/2. The amplitude Z p is defined in 
terms of the point-to-point pseudoscalar propa- 
gator through 

(0|P(t)P(0)t|0) 2Z P coshm polo (i - L t /2).(4) 

In Fig. H we plot a typical example of the ratio 
of propagators in (||). Data for the ratio for the 
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Figure 3. Effective mass for heavy-light meson. 
Parameters and meaning of symbols are the same 
as in Fig. 0. 



wall source exhibits a smooth behavior with small 
errors of 1% or less for any value of the hopping 
parameter we took. Thus the constant Z A l p can 
be reliably determined. 

For an extraction of the amplitude Z p we em- 
ployed a single hyperbolic cosine fit to the point- 
to-point pseudoscalar propagator at the time of 
the Symposium. We observed that this is prob- 
lematical when the heavy quark hopping param- 
eter becomes small because of a rapid increase of 
fluctuations. This is illustrated in Fig. || which 
shows the effective mass of the heavy-light meson 
at our smallest value of Kq: no clear signal of 
plateau can be seen for the point-to-point propa- 
gator (open circles). 

An alternative procedure, facilitated by the 
fact that signals are much better for the wall 
source, is to apply a simultaneous fit to the point- 
to-point and point-to-wall propagators assuming 
a common mass. We find that the fit becomes 
much more stable and that the fitted value of 
Z p over subsets of configurations fluctuates much 
less. We therefore adopt this procedure to extract 
the value of Z p . An example of the fitted value 
and the error for mass is shown in Fig. |^ (solid 
and dashed lines). 

At = 6.3 we may compare our results for 
lattice value of fp^/mp with those of Ref. 0. 
Applying the same analysis procedure as ours to 
their data we find that the two sets of results are 
consistent within the error. 
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Figure 4. f p as a function of inverse kinetic 
mass 1/mp in physical units. 

We employ the tadpole-improved one-loop 
renormalization factor Z4, already encountered 
for f„, to convert lattice results for fpyjmp to 
those in the continuum. Since this factor refers 
to the case of zero quark mass, an uncertainty of 
0(a s rriQa) remains in our final result. In Fig. || 
we plot the scale invariant quantity 



a s (mB) J 



(5) 



as a function of the inverse of meson mass 1/mp 
for both /3=6.1 and 6.3, where we take the mea- 
sured kinetic mass for mp. We also plot the static 
values taken from Refs. ^J,[?J for comparison. Con- 
version to physical units is made with the lattice 
spacing listed in Table |. 

We observe that the results at the two values 
of j3 agree with each other, and they seem to con- 
tinue smoothly to the static value. Errors in the 
kinetic mass and statistical fluctuations arising 
from point-to-point propagators in the values for 
$p, however, are still large to tell conclusively 
whether large 0(mQa) effects are absent for the 
heavy- light meson decay constant. We remark 
that this result is obtained when $p is plotted 
with respect to the measured kinetic mass. 

We fit results for $p toa quadratic form $p = 
co(l + ci/mp + 02 /mp) to obtain the values at 
the physical B and D mesons. Our results for fp, 
Jb b ) /d and fu. are given in Table ||. Statistical 
errors include uncertainties in the kinetic mass to 
the extent that a jackknife analysis eliminating 



Table 2 

Decay constant and B parameter for B and D 
mesons. First error is statistical and second that 



of scale determined from m p . 




6.1 


6.3 


1b (MeV) 


192(5)(10) 


182(10)(12) 


fB, (MeV) 


228(4)(12) 


215(4)(14) 


/b, / !b 


1.19(3) 


1.18(6) 


Id (MeV) 


206(5)(11) 


216(9)(14) 


Id. (MeV) 


237(4)(13) 


240(4) (16) 


Jd 3 1 Id 


1.16(3) 


1.11(5) 


B b 


0.895(47) 


0.840(60) 


Bb 3 


0.889(24) 


0.878(32) 



single gauge configuration at a time is applied to 
the fitting of $ P . 

6. B meson B parameter 

We also calculate the B meson B parameter 
from a ratio of three- and two-point functions and 
the AB = 2 weak operator fixed at the origin of 
lattice. One-loop renormalization factors evalu- 
ated with the tadpole-improved coupling ay (1/a) 
are used for relating the lattice value to that in 
the continuum. Our results for Bp and Bp s at 
/x=5GeV are tabulated in Table ||, which are con- 
sistent with results reported previously 
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